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This  thesis  describes  the  use  of  the  NASTRAN  program 
(hovel  17.0)  in  the  calculation  of  flexible  wing  airloads 
and  stresses.  The  problems  of  interfacing  the  aerodynamic 
and  structural  models  are  discussed  with  assumptions  need¬ 
ed  to  solve  them.  Two  different  methods  of  transferring 
the  aerodynamic  forces  into  a  structural  load  vector  are 
presented . 

A  Direct  Matrix  Abstraction  Program  was  developed  for 
use  in  these  calculations.  This  sequence  is  limited  to 
the  use  of  one  doublet- lattice  panel  for  the  .aerodynamic 
model  and  requires  several  new  data  blocks  for  execution. 
Although  these  data  blocks  were  input  with  the  bulk  data 
deck  to  test  the  sequence,  a  new  preliminary  module  has 
been  constructed  to  build  them  internally. 

The  results  of  this  investigation  show  that  NASTRAN 
can  be  used  in  calculating  the  airloads  and  stresses  for 
a  flexible  wing.  This  new  sequence  has  extended  the  capa¬ 
bility  of  NASTRAN  to  allow  for  the  application  of  inter¬ 
nally  generated  steady  airloads  to  the  structural  model. 
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CALCULATION  OF  AIRLOADS  FOR  A  FLEXIBLE 
WING  VIA  NASTRAN 

I  Introduction 

The  calculation  of  flexible  aircraft  airloads  and  struc¬ 
tural  stresses  requires  the  interfacing  of  both  structural 
and  aerodynamic  theories.  The  NASTRAN  structural  analysis 
program  includes  both  of  these  types  of  theories.  Tne 
purpose  of  this  study  was  to  demonstrate  the  feasibility  of 
using  this  program  for  these  calculations. 

The  NASTRAN  program  uses  a  finite  element  structural 
model.  A  wide  variety  of  element  types  are  available,  thus 
allowing  the  construction  of  simple  to  very  complex  models. 
The  aerodynamic  theory  chosen  for  this  study  uses  doublet- 
lattice  panels  to  model  the  lifting  surfaces.  Other  aero¬ 
dynamic  theories  are  available  for  use  in  the  NASTRAN  pro¬ 
gram  but  were  not  considered. 

A  new  instruction  seguence  Known  as  a  Direct  Matrix 
Abstraction  Program  (DMAP)  was  developed  for  use  within 
NASTRAN  to  calculate  flexible  aircraft  airloads  and  struc¬ 
tural  stresses.  This  new  DMAP  5-  'quonco  is  limited  to  a 
structural  model  consisting  of  only  a  wing.  The  correspond¬ 
ing  aerodynamic  model  is  limited  to  a  single  aerodynamic 
panel  to  model  the  lifting  surface.  This  new  seguence 
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extends  the  capability  of  the  NAS TRAN  program  by  creating 
the  ability  to  apply  the  internally  generated  steady  air¬ 
loads  to  the  structural  model.  The  potential  to  do  this 
has  existed  since  the  inclusion  of  aerodynamic  theories  in 
the  NASTKAN  program. 

The  calculation  of  flexible  wing  airloads  by  this  DMAP 
seguence  is  an  iterative  process.  The  initial  airloads  are 
calculated  based  upon  a  rigid  wing  shape  at  an  nngle-of- 
attack.  These  airloads  are  applied  to  the  structural  model 
and  the  deflections  calculated.  The  deflected  structure 
causes  the  airloads  to  change,  which  will  in  turn  change 
the  deflection  of  the  structure.  This  cyclic  process  is 
continued  until  the  incremental  changes  in  the  airloads  and 
the  displacements  are  negligible. 

The  wing  model  chosen  to  demonstrate  the  new  DMAP  se¬ 
quence  was  first  tested  for  solution  convergence  of  the 
cyclic  process.  The  test  was  performed  using  the  rigid 
format  for  static  analysis  interfaced  with  the  Woodward 
aerodynamic  program.  The  results  of  this  test  were  also 
used  as  an  alternative  solution  to  check  the  reasonableness 
of  the  solution  generated  by  the  new  DMAP  sequence. 

The  potential  users  for  this  new  DMAP  seguence  are  the 
Structures  Branch  of  the  Fliqht  Dynamics  Lab  ( AFr DL/FB )  and 
San  Antonio  Air  Logistic  Center  (SAALC).  Tne  SA  ,LC  has  had 
a  NASTRAN  structural  model  of  half  of  a  T-37  aircraft  since 
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1976.  This  model  was  in herded  to  be  used  to  calculate 
stresses  developed  by  an  aerodynamic  loading  of  the  model . 

A  workable  automated  method  of  generating  airloads  for  this 
model  has  not  been  found.  The  AFFDL/FB  has  developed  two 
programs  to  generate  loads  for  a  wing  structural  model 
from  airloads  calculated  by  an  aerodynamic  program.  Both 
organizations  can  use  the  new  DMAP  sequence. 


3 


II  Background 


| 

)  A  current  computer  program  for  calculating  the  loads 

t 

! 

;  on  a  flexible  aircraft  is  Flexloads.  This  program  includes 

several  theories  to  calculate  aerodynamic  forces  for  veloci¬ 
ties  ranging  from  subsonic  to  supersonic.  These  aerodynamic 
theories  are  the  Woodward  Finite  Element  Subsonic-Supersonic 
Steady  State  Theory,  Doublet-Lattice  Steady  and/or  Unsteady 
Theory,  Steady  and/or  Unsteady  Modified  Newtonian  Theory, 
and  finally  Subsonic  and  Supersonic  Kernal  Function  Theory. 
The  structures  are  modeled  in  this  program  as  slender  beams 
or  as  built-up  plates.  As  a  result.  Flexloads  cannot  use 
all  the  different  elements  available  with  NASTRAN. 

The  AFFDL/FB  has  developed  two  programs  to  generate 
equivalent  loads  for  a  NASTRAN  wing  model  from  airloads 
calculated  by  an  aerodynamic  program.  These  two  programs 
are  SPLOADS  and  BEAMPROG.  These  programs  handle  the  inter¬ 
facing  problem  of  applying  airloads  to  a  structural  model. 
These  programs  could  be  used  as  part  of  the  iterative  pro¬ 
cess  necessary  in  the  calculation  of  flexible  aircraft 
j  airloads.  However,  another  program  would  have  to  be  devel- 

t 

oped  to  introduce  structural  deformation  effects  into  the 

j 

calculation  of  airloads  before  this  process  could  be  em¬ 
ployed  . 
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Ill  Inter!  acl  aq  Problems 

The  problems  of  interfacing  an  aerodynamic  wing  model 
with  its  corresponding  structural  model  are  discussed  in 
this  section.  The  two  models  solve  two  different  problems, 
so  the  mesh  sizes  used  for  the  two  models  will  be  different, 
which  causes  many  of  the  interfacing  problems.  The  inter¬ 
facing  is  required  for  the  application  of  airloads  to  a 
■:l  riii'lur.il  model  and  for  I  he  ini  roduel  ion  of  si  ruclural 
deformation  into  the  calculation  of  airloads.  In  the  fol¬ 
lowing  these  problems  are  discussed  in  general  and  as  they 
specifically  apply  to  the  new  DMAP  sequence. 

The  structural  models  of  most  wings  only  include  the 
main  structural  elements,  the  wing  box.  This  means  the 
planform  of  the  aerodynamic  model  is  larger  than  that  of  the 
structural  model.  Because  of  this  difference,  the  aerodyna¬ 
mic  model  will  generate  loads  for  which  the  corresponding 
structural  elements  do  not  exist.  These  airloads  have  to 
be  handled  as  special  cases  for  translation  and  application 
to  the  structural  model.  The  difference  in  planform  sizes 
causes  problems  in  the  calculation  of  new  airloads  based 
on  structural  del ormation .  The  structural  deformation  changes 
the  local  slope  of  an  aerodynamic  box  which  is  used  in  the 
calculation  of  the  down  wash.  However,  the  deformed  structure 
necessary  for  this  calculation  does  not  exist  in  areas  for 


the  aerodynamic  model  lying  outside  the  wing  tox.  A  simple 
solution  to  this  difference  in  planform  sixes  Is  to  force  the 
two  to  be  the  same.  As  a  result  the  structural  model  must 
include  the  secondary  wing  structure.  Since  this  secondary 
wing  structure  only  transfers  the  loads  to  the  wing  box,  this 
structure  does  not  have  to  be  modeled  accurately  except  in 
size  of  the  planform. 

Even  when  these  planforms  are  the  same,  the  aerodynamic 
I  <>  i  ;  have  fo  bo  I  r.msfoniioj  inlo  a  )oa<J  vool  nr  for  I  lie  struc¬ 
tural  model.  The  doublet-lattice  method  used  in  NASTRAN  gen¬ 
erates  a  delta  pressure  coefficient  AC  which  is  the  difference 

P 

between  the  pressure  on  the  upper  and  lower  surfaces  at  the 
aerodynamic  points  for  each  box.  For  simplicity  the  new  DMAP 
sequence  applies  the  load  to  the  upper  surface  of  the  struc¬ 
tural  model.  Assumptions  on  the  AC^  distribution  are  required 
to  transform  the  aerodynamic  forces  into  a  load  vector.  Two 
methods  based  upon  a  different  set  of  assumptions  for  the 
distribution  have  been  implemented. 

In  the  first  method  the  AC  value  at  the  structural  ele- 

P 

ment  centroid  is  found  by  a  surface  spline  interpolation  and 

then  assumed  constant  over  the  element.  Because  AC  is  con- 

P 

stant,  multiplication  by  the  dynamic  pressure  and  the  element 
area  will  yield  the  total  load.  For  basicaJIy  parallelogram 
shaped  elements,  one  quarter  of  this  total  load  .s  now  distri¬ 
buted  to  each  of  the  element's  nodes.  This  is  done  for  each 

membrane  element  on  the  upper  surface,  result  inq  in  a  global 

0 


expression  for  the  load  vector. 


In  ';he  second  transformation  method  more  variation  in 
ACn  distribution  is  accounted  for.  After  finding  interpolat- 
ed  corner  values  for  each  element,  a  bilinear  variation  is 
assumed  for  the  distribution  over  two  triangular  subregions. 
See  Figure  1. 


Figure  1.  Division  of  the  Upper  Surface  Quadralateral 
Membrane  Elements  into  Triangular  Elements. 

The  bilinear  pressure  over  each  of  the  triangular  subelements 
is  assumed  to  £>a 

P  (X,  Y)  =  q  Li  ACp. 

where  q  is  the  dynamic  pressure,  Li  are  the  shape  functions 


for  a  three  noded  trianqular  element,  and  ACn  is  a  vector 

Hi 

of  nodal  values  of  ACp.  The  equivalent  nodal  forces  are 
calculated  as  usual  by: 

Fi  -  4  PLidA  =  J'a  dLj  CPiLidA 

When  integrated  this  becomes 


Those  subolomont  expressions  are  now  assembled  into  a  global 
expression  for  the  load  vector. 

The  doublet- latt ice  method  calculates  the  aerodynamic 
forces  normal  to  the  panel .  When  the  wing  is  at  an  angle- 
of-attack  the  orientation  of  the  lift  vector  and  the  normal 
force  vector  will  differ  by  this  angle.  The  new  DMAP  sequence 
uses  the  normal  force  instead  of  lift  in  calculating  the  co¬ 
efficients  of  lift,  pitching  moment  and  rolling  moment.  The 
resulting  errors  in  these  coefficients  will  be  minimal,  if 
small  angles-of-attacks  are  used. 
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TV  Procedures 

An  intermediate  complexity  wing  structural  model  was 
obtained  from  AFFDL/FB  to  demonstrate  the  new  DMAP  sequence. 

A  description  of  the  model  used  is  contained  in  Appendix  B. 
The  wing  model  was  first  tested  by  using  the  Woodward  pro¬ 
gram  USSAERO  to  generate  the  airloads  used  in  the  iterative 
process.  This  test  checked  the  model  for  solution  conver¬ 
gence  nnd  provided  an  nUernnt  i  ve  solution.  The  new  DMAP 
sequence  was  then  used  to  resolve  the  problem. 

This  section  covers  first  the  interfaced  NASTRAN-USSAERO 
sequence  and  next  the  new  DMAP  sequence.  The  nev  sequence 
is  explained  relative  to  functions  it  performs,  new  data 
blocks  it  requires,  and  in  how  it  would  use  the  new  module. 

interfaced  NASTRAN-USSAERO  Sequence 

The  static  analysis  rigid  format  of  NAS TRAN  and  USSAERO 
are  interfaced  in  a  manual  execution  of  the  iterative  pro¬ 
cess.  The  output  from  USSAERO  is  used  to  create  constant 
pressure  loads  on  upper  surface  structural  elements  and  then 
the  displacements  from  NASTRAN  are  used  to  change  the  camber 
surface  for  USSAERO  in  this  process. 

The  doublet-lattice  method  in  NASTRAN  is  a  lifting  sur¬ 
face  theory,  therefore  USSAERO  which  includes  thickness  ef¬ 
fects  approximates  the  thin  sheet  by  using  the  planar  option 
and  the  smallest  allowable  winq  thickness. 

<1 


USSAERO  generates 


values  for  pressure  coefficient  Cp  at  the  center  of  each 
aerodynamic  box  for  both  the  upper  and  lower  surfaces.  A 
further  explanation  of  USSAERO  is  contained  in  Reference 

8. 

The  mesh  size  for  the  aerodynamic  model  was  chosen 
to  avoid  interface  problems  by  exactly  matching  that  for 
the  structural  model.  By  this  matching  the  need  to  trans¬ 
form  the  airloads  for  structural  application  is  avoided. 

The  l  dI  a  L  demon l  a. L  Load  now  lel.iLe:;  lo  l  he  ae  l  ody  liain  i  c 
force  produced  hy  the  corresponding  box. 

The  boundary  conditions  for  the  static  analysis  were 
chosen  to  take  out  the  rigid  angle-of-attack  displacements, 
leaving  only  total  deformation  in  the  looping  process.  To 
accomplish  this,  all  degre es-of -f reedom  for  the  wing  root 
grid  points  were  fixed.  For  each  loop  pass,  the  new  total 
deformation  were  calculated  from  the  changing  total  airloads 

Two  conversion  programs  were  reguired  to  interface 
NASTRAN  and  USSAERO  for  the  iterative  process.  Because  of 
the  matched  meshes,  these  essentially  just  changed  the 
output  from  one  into  the  proper  format  for  input  into  the 
other.  The  first  program  converts  the  AC  to  pressure  loads 
To  do  this  the  delta  coefficient  of  pressure  for  each  box 
is  found  and  then  multiplied  by  the  dynamic  pressure  g. 

These  pressure  loads  are  used  to  create  PL0AD2  c..rds  for 
use  in  the  NASTRAN  bulk  data  deck.  The  second  p.  ograin 
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|  converts  the  total  deformation  to  a  specification  of  a  new 

i. 

t 

f  camber  surface  fot  the  aer^dynanic  model.  The  structural 


displacements  are  averaqed  to  find  a  midsurface,  which  also 


defines  the  camber  surface  used  by  USSAERO.  The  position 


of  these  two  programs  in  the  iterative  process  is  shown  in 


Figure  2. 


Figure  2.  INTERFACED  NASTRAN-USSAERO  Sequence 


This  manual  execution  of  the  iterative  process  is  in¬ 
efficient  for  the  following  reasons.  The  equations  have  to 
be  reformed  with  every  loop  pass  when  only  the  right  hand 
side  of  the  equations  change.  Next,  the  process  generates 
large  amounts  of  intermediate  data  that  is  not  used.  Also, 
the  precision  of  the  intermediate  solutions  are  decreased, 
because  the  data  passed  from  one  program  to  the  next  must 
be  truncated  to  specific  formated  field  sizes.  Finally, 
the  process  requires  sequential  jobs  to  be  run  with  each 
pass,  requiring  excessive  amounts  of  user  time. 

The  DMAP  Sequence 

The  new  DMAP  sequence  calculates  the  airloads  and 
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structural  stresses  for  a  flexible  winq  model  .  ,phe  airloads 
used  in  this  sequence  are  calculated  by  the  doublet-lattice 
method  in  NASTRAN.  Ry  us i nq  a  reduced  frequency  of  approxi¬ 
mately  zero,  this  unsteady  method  reduces  to  a  steady  vortex- 
lattice  method.  Values  of  AO  are  produced  for  the  center 
of  the  quarter-chord  for  each  box.  A  further  explanation 
of  this  aerodynamic  method  is  found  in  References  2,  3  and 
G. 

The  new  DMAP  sequence  is  explained  by  major  t unctions 
in  the  following  section.  The  total  capability  lias  previous¬ 
ly  existed  for  NASTRAN  to  calculate  the  flexible  wing  airloads, 
except  for  the  ability  to  transform  the  steady  airloads  into 
a  structural  load  vector.  io  test  the  new  DMAP  sequence, 
the  data  blocks  necessary  to  build  this  transformation  matrix 
were  input  directly  with  the  bulk  data  deck.  A  discussion 
of  these  new  data  blocks  follows  the  explanation  of  the  se¬ 
quence.  A  new  module  was  designed  that  will  build  this  trans¬ 
formation  matrix  internally.  The  functions  of  the  new  module 
are  explained  in  the  section  following  the  data  block  descrip¬ 
tions  . 

Several  restrictions  are  required  on  models  used  with 
the  new  DMAP  sequence  as  it  presently  exists.  The  plan- 
forms  for  the  aerodynamic  and  structural  models  are  required 
to  be  the  same.  This  is  necessary,  as  discussed  in  Section 
III,  so  that  the  airloads  can  be  applied  to  correspondinq 

structural  elements  without  consideration  of  special  cases. 
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Next  the  aerodynamic  model  is  restricted  to  one  doublet- 
lattice  panel  to  simplify  bookkeeping.  Finally,  an  analy¬ 
sis  set  must  be  specified  that  contains  only  the  Z-direction 
degree-of-f reedom  for  the  wing  grid  points  exclusive  of  the 
wing  root.  This  limitation  is  reguired  for  correct  calcu¬ 
lation  of  the  airloads. 

Functions  of  the  New  DMAP  Sequence .  The  new  DMAP  se¬ 
quence  will  be  described  here  in  general  terms  by  the  major 
I  unci,  ions  it  performs.  'I' lie  complete  .Listing  lot  l  lie  UMA1’ 
sequence  is  contained  in  Appendix  A.  The  sequence  consists 
largely  of  a  combination  of  rigid  formats  used  for  static 
analysis  and  flutter  analysis.  The  sequence  accomodates 
any  of  the  structural  element  types  available  for  use  with 
the  static  analysis  rigid  format,  but  does  not  allow  the 
use  of  the  down  wash  factors  due  to  extra  points  available 
in  the  flutter  analysis  rigid  format. 

The  main  functions  of  the  new  DMAP  sequence  will  bo 
presented  in  sequential  order.  The  first  portion  of  the 
sequence  sets  up  the  tables  necessary  to  describe  the  struc¬ 
tural  model,  including  the  element  connection  table  (ECT)  and 
the  grid  point  list  (GPL) .  These  tables  are  used  in  the  as¬ 
sembly  of  the  unconstrained  global  stiffness  and  mass  matrices 
for  the  system. 

The  constraints  for  subcase  one  are  now  applied  to  these 
global  matrices.  These  constraints  are  constructed  so  as 
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to  leave  only  a  rigid  body  pitch  mode.  The  coast  rained  ma¬ 
trices  are  then  partitioned  down  to  the  analysis  set  speci¬ 
fied  by  the  ASET  card. 

A  real  eigenvalue  analysis  is  performed  on  the  analysis 
set  to  calculate  the  displacements  due  to  this  rigid  pitch 
mode.  The  mass  matrix  is  only  used  in  this  sequence  for 
the  eigenvalue  analysis,  therefore  a  lumped  mass  matrix  is 
sufficiently  accurate.  This  riqid  pitch  mode  shape  is  now 
mill  I  iplied  by  an  i.npul  parameter  AOAP  to  speeiiy  displace¬ 
ments  for  an  initial  angle-of-uttack.  This  procedure  was 
necessary  because  the  NASTU’AN  program  does  not  have  t  hi' 
capability  to  set  the  model  at  an  initial  angle-of-attack. 

The  camber  is  now  added  to  the  displacements  due  to  angle- 
of-attack. 

The  next  function  is  the  formation  of  the  transforma¬ 
tion  matrix  GTKA  which  transforms  the  displacements  to  values 
of  slope  at  the  aerodynamic  points.  These  slopes  .are  used 
to  specify  the  strengths  for  the  steady  down  wash  vector 
D1JK .  The  aerodynamic  influence  matrix  AJJL  for  the  doublet- 
lattice  method  is  calculated  next.  The  inverse  of  this  ma¬ 
trix  when  multiplied  by  the  down  wash  vector,  yields  the 
values  for  AC^  at  each  of  the  aerodynamic  points.  These 
ACp  values  are  for  the  total  pressure  differential  across 
♦■he  wing. 

The  next  functions  of  the  sequence  are  those  the  new 


moduli'  will  perform  when  incorporated  inLo  NASTNAN.  1'Jie 

L  runs  format  ion  matrix  TSTA  is  :-a J  cula  l  (’d  to  transform  t  ho 

vector  of  AC  values  to  an  unconstrained  global  load  vector 
P 

PG.  This  load  vector  will  only  contain  non-zero  forces  for 
the  Z-direction  deqree-of- f reedom  for  the  upper  surface 
grid  points.  The  matrices  necessary  for  the  calculation 
of  the  coefficients  of  lift,  pitching  moment,  and  rolling 
moment  are  also  built  in  this  section  of  the  sequence. 

The  hi  >i  in  da  r  y  i’oikJ  i  I  i  on::  are  now  el  ian<  |e<  j  I  i  >  I  ho;  ;e  spec¬ 
ified  by  subcase  two,  which  constrains  all  the  deqrees-of- 
f reedom  for  the  winq  root  qrid  points.  The  matrix  data 
blocks  necessary  for  a  static  deflection  calculation  are 
setup,  with  this  boundary  condition.  The  displ  jcemonts 
calculated  will  be  from  the  midplane  of  the  winq  pitched 
at  an  anglo-of-attack .  Those  matrices  are  partitioned  down 
to  the  same  analysis  set  used  to  calculate  the  airloads, 
hy  doing  this,  the  transformation  matrices  TSTA  and  GTKA 
previously  calculated  can  be  used  in  the  iterative  loop 
to  calculate  the  flexible  wing  airloads. 

The  next  function  of  the  sequence  is  the  loop  to  cal¬ 
culate  the  flexible  wing's  displacements,  static  load  vector 
and  airloads  by  the  iterative  process.  This  loop  calculates 
an  incremental  displacement  vector  due  to  the  initial  load 
vector.  This  incremental  displacement  vector  is  then  used 
to  calculate  an  incremental  airload  vector  which  is  transformed 
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lo  an  incremental  load  vector.  The  iterative  process  is 
continued  until  the  norm  of  the  incremental  displ acoment 
vector  is  less  than  a  prescribed  input  parameter  (thooP). 

The  incremental  displacements,  incremental  load  vectors 
and  incremental  airloads  are  summed  with  each  pass  through 
the  loop.  These  sums  give  the  final  displacement  of  the 
analysis  set,  the  final  load  vector  and  the  final  airloads 
for  the  flexible  wing. 

The  nexl  purl  of  l  h(>  sequence  i(  '<  •<  >ve  r: ;  I  1 1<  >  ( •<  imp  I  < '  l  ( ' 
solution  for  the  analysis  sot  and  the  omitted  deqrees-of- 
freedom.  iTom  tne  complete  solution,  the  stresses  are  now 
calculated  for  any  elements  specified  in  subcase  two  of  the 
case  control  deck. 

The  coefficients  of  lift,,  pitching  moment  and  rolling 
momenL  due  to  both  the  final  airloads  and  the  final  load 
vector  are  calculated  next.  These  coefficients  are  used  as 
a  check  on  errors  induced  when  the  airloads  are  transformed 
to  a  load  vector.  The  final  section  of  the  sequence  contains 
error  messages  called  when  a  specific  error  is  found  during 
the  execution  of  the  sequence. 

New  Data  Blocks  and  Parameters .  The  data  blocks  re¬ 
quired  by  the  new  DMAP  sequence  wore  created  and  input  direct¬ 
ly  with  the  bulk  data  deck.  All  of  the  input  matrices  except 
the  camber  vector  will  be  created  by  the  new  moduli'  when  it 
is  incorporated  into  NASTRAN.  Some  extra  steps  were  required 
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in  Lhu  sequence  to  process  Lliese  input  mu  trices  unci  will  be 
removed  when  the  module  is  added. 

The  new  parameters  that  ur?  needed  for  this  sequence 
are  listed  with  a  brief  explanation  in  Table  I.  A  list  of 
the  new  input  matrices  for  the  sequence  is  given  in  Table 
II.  The  major  function  of  these  matrices  is  the  transfor¬ 
mation  of  the  aerodynamic  forces  to  the  structural  load 
vector  PC. 

Thin  loud  vcf’lor  cun  bo  ore.ilod  by  one  of  I  he  two 
methods  discussed  in  Section  III.  The  first  method  gener¬ 
ates  the  transformation  matrix  by 

[TSTA]  =  q  [TTTTT]  [SASEL]  [GTaKT]T 

for  use  in 

|PG}  =  [TSTA]  |QKKK} 

which  PC  is  the  unconstrained  load  vector  and  QKKK  is  the 
vector  of  ACp  values  at  the  aerodynamic  points.  The  trans¬ 
formation  matrix  GTAK  for  this  method  interpolates  the  QKKK 

matrix  to  a  vector  of  AC  values  at  the  structural  element 

P 

centroids . 

The  second  method  generates  the  transformation  matrix 
TSTA  in  a  slight  different  way.  The  matrix  is  now 
[TSTA]  =  q  [TTSAS]  [GTAKT]T 

where  the  matrix  TTSAS  is  the  assembled  global  interpolation 
matrix  for  a  consistent  load  vector  for  triangular  elements 
and  the  matrix  GTAKT  interpolates  the  same  QKKK  matrix  used 
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Table  I. 


List  of  Now  Parameters. 


NAML 


DESCRIPTIONS  OP  PARAMETER 


AOAP 


Angle  of  attack  parameter 
angle  of  attack;  Complex 


,  used  to  set 
( AOAP ,0.0) 


QC  Dynamic  pressure,  used  for  calculation 

of  loads;  Complex  (Q,0.0) 


u 

Dynamic  pressure,  used  in  calculation 
of  aerodynamic  coefficients;  Real 

X 

n 

Reference  chord  length,  used  in  calcu¬ 
lation  of  aerodynamic  coefficients; 

Real 

WAREA  Plan form  area  of  the  structural  model; 

Real 


ICAMB 


Flag  for  camber  vector,  only  input  if 
camber  to  be  specified;  Integer >0 


ELOOP 


Used  to  specify  convergence  criteria 
for  the  norm  of  the  incremental  dis¬ 
placements  Real  >0  or  as  a  flag  to 
specify  only  rigid  ACp  calculations 
Real  <0 
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Table  II.  List  of  New  Matrices 


NAME 

DESCRIPTIONS  OF  MATRICES 

OTAKT 

Surface  spline  interpolation  matrix  for  trans¬ 
formation  of  ACp's  at  the  aero  points  to  ACp's 
at  structural  points;  real,  general 

SASEL 

Contains  areas  of the  surface  membranes  to  which 
loads  are  applied;  real,  diagonal 

TT'PTT 

Transformal'  i  on  matrix  to  distribute  1/d  of  the 
load  on  Llie  membranes  Lo  Lheir  nodes;  real, 
general 

TVS  AS 

Assembled  global  interpolation  matrix  for  con¬ 
sistent  load  vector,  triangular  elements,  re¬ 
places  SASEL  and  TTTTT  lor  linear  pressure- 
distribution;  real,  general 

CAMBER 

Contains  camber  surface  displacements  i or  the 
analysis  set  grid  points;  real,  vector 

SAS 

A  matrix  containing  one's  where  areas  are  in 

SKJ  matrix,  imaginary  part  zero;  complex,  gen¬ 
eral 

Unit  vector  of  length  equal  to  number  of  aero¬ 
dynamic  boxes;  real,  vector 

_  _  . 

VGLS 

A  vector  with  ones  in  the  proper  positions,  so 
that  when  multiplied  by  PG  it  sums  the  loads; 
real,  vector 

VGMS 

A  matrix  with  ones  in  the  proper  positions,  so 
that  when  multiplied  by  PG  it  creaLes  a  vector 
of  the  non-zero  loads;  real,  general 

VTMA 

Contains  the  moment  arms  for  the  aerodynamic 
pitching  momer t  ibout  the  X-axis;  real,  vector 
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Table  II 


List  of  New  Matrices  (Cont.' 


NAME 

DESCRIPTIONS  OF  MATRICES 

VTRA 

Contains  the  moment  arms  for  the  aerodynamic 
rollincr  moment  about  the  Y-axis;  real,  vector 

VTM 

Contains  the  moment  arms  for  the  pitching 
moment  due  to  the  transformed  loads  about  the  . 
Y-axis;  real,  vector 

VTR 

Contains  the  moment  arms  for  the  rolling 
moment  due  to  the  transformed  loads  about  the 
X-axis;  real,  vector 
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in  the  first  method  to  a  vec  :or  of  nodal  AC  values.  The 

P 

equation  for  the  unconstrained  consistent  load  vector  is 
jPG}  =  [TSTA"  jQKfK } 

where  PG  and  QKKK  are  defined  the  same  as  in  the  first 
method . 

formation  of  the  interpolat ion  matrix  OTAKT  is  the 
same  for  both  methods.  The  NASTRAN  subroutines  source 
decks  used  in  the  surface  spline  interpolation  were  used 
l.o  build  this  matrix  in  a  temporary  external  proyram.  By 
using  the  same  subroutines  the  transposed  interpolation 
matrix  will  be  the  same  as  the  one  generated  by  the  module 
when  incorporated.  Since  the  number  of  significant  figures 
is  limited  to  the  size  of  the  data  fields  used  fo\  input, 
the  large  field  bulk  data  cards  should  be  used. 

The  camber  vector  listed  on  Table  11  is  used  to  speci¬ 
fy  the  camber  surface  of  the  wing,  when  required.  This 
vector  contains  the  Z-direction  camber  displacements  for 
the  upper  surface  grid  points  in  the  analysis  set.  Zeros 
are  entered  for  the  lower  surface  grid  points.  In  addition, 
if  the  wing  is  twisted  this  is  added  into  the  camber  vector 


in  a  similar  manner 


Module  Description.  The  new  module  when  incorporated 
into  the  NASTRAN  program  will  be  designated  by  the  name  G  r  AS . 
This  module  was  not  incorporated  into  the  NASTRAN  program 
because  the  link-editor  routine  for  the  program  was  not  func¬ 
tional.  This  module  will  build  internally  all  of  the  matrices 
listed  in  Table  II,  except  for  the  camber  vector.  The  module 
will  be  called  in  the  new  DMAP  sequence  prior  to  the  calcu¬ 
lation  of  the  initial  AC  values.  The  exact  posit  Lon  of 

r 

this  new  module  and  Liie  step:;  it  will  replace  in  the  sequence 
are  shown  in  Appendix  A.  The  fortran  source  program  for  this 
module  is  contained  in  Appendix  C.  This  source  program  has 
not  been  tested  wish  NASTRAN,  but  it  should  function  correct¬ 
ly  when  incorporated. 

The  following  restrictions  are  required  when  using  this 
module.  Since  only  the  upper  surfaces  is  loaded,  the  struc¬ 
tural  model  is  restricted  to  CQMEM2  membrane  elements  on  that 
surface.  Also,  the  aerodynamic  model  is  restricted  to  one 

doublet-lattice  panel  which  must  lay  in  the  X-Y  plane. 

\ 

The  DMAP  calling  sequence  for  this  module  will  be 

GIAS  ACPT ,TOSEL,ECTA ,GPLA ,BGPA/SAS , 

GT AK , TTTTT , S  AS  E  L , VGLS , VGMS , VTMA , 

VTRA , VTM , VTR , AI DMIT/$ 

The  five  data  blocks  before  the  first  slash  mark  are  the 
required  input  tables  for  this  module.  These  input  tables 
are  listed  in  Table  III  with  a  short  descr ipt ron .  The  e- 
leven  output  data  blocks  listed  between  the  slash  marks  are 
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the  same  matrices  described  in  Table  II. 


Table  III.  Input  Tables  for  the  Module 


NAME 

DESCRIPTION 

ACPI 

Aerodynamic  connection  property  table  con¬ 
tains  the  value  for  the  total  number  of 
aero  boxes . 

TOSEL 

Sorted  list  of  CQMEM2  elements  to  which 
loads  will  bo  aoDlied.  in  uscendina  order. 

KOTA 

Element  connection  table  for  both  struc¬ 
tural  o.lemenls  and  aero  boxes. 

OPT,  A 

Grid  point  numbers  listed  in  order  of 
appearance  in  the  BGPA  table. 

BGPA 

List  of  coordinates  for  both  structural 
arid  points  and  arid  points  of  aero  boxes. 

The  TOSEL  table  is  specifically  required  by  the  new 
module  and  must  be  input  with  the  bulk  data  deck.  This 
table  lists  the  CQMEM2  element  identification  numbers  for 
the  elements  on  the  upper  surface.  These  are  the  elements 
the  aerodynamic  force  will  be  applied  to  when  transformed 
to  a  load  vector.  These  numbers  must  be  listed  in  order 
from  lowest  to  highest. 

This  module  builds  the  data  blocks  necessary  for  load 
application  by  method  one.  For  this  method,  the  interpola¬ 
tion  matrix  GTAK  is  built  by  a  surface  spline  through  the 
aerodynamic  points  located  at  the  center  of  the  quarter  chord 
of  each  box  and  the  centroids  of  the  membrane  elements.  This 


module  also  builds  the  data  blocks  necessary  for  the  calcu¬ 
lation  of  the  aerodynamic  coefficients.  All  of  these  data 
blocks  are  built  with  the  origin  as  the  reference  point  for 
calculating  the  coefficients. 
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V  Results 

The?  new  DMAP  sequence  was  successful  in  calculating  the 
airloads  and  stresses  for  t  flexible  wing.  Several  sample 
cases  were  run  and  are  presented  in  this  section.  Results 
were  obtained  for  two  different  aerodynamic  models  and  for 
two  methods  of  airload  transformation.  These  are  compared 
with  the  solution  obtained  by  the  NASTRAN-USSAERO  sequence. 
Comparisons  of  rigid  wing  AC^  distribution,  errors  caused 
by  the  assumption  used  for  load  transformation,  wing  tip 
deflections,  and  the  convergence  of  the  iterative  processes 
are  presented.  A  summary  of  the  cases  compared  are  contain¬ 
ed  in  Table  IV. 

The  doublet-lattice  panel  is  an  unsteady  thin  airfoil 

aerodynamic  model,  that  for  zero  frequency  reduces  to  a 

steady  vortex  lattice  panel.  The  new  DMAP  sequence  only 

uses  the  steady  part  of  this  panel,  so  this  part  was  checked 

against  USSAERO  for  accuracy.  The  chordwise  AC  values  cal- 

P 

culated  for  a  rigid  wing  are  shown  in  Figure  3.  The  steady 
part  calculates  lower  values  of  ACp  near  the  leading  edge 
and  slightly  higher  near  the  trailing  edge.  The  AC^  dis¬ 
tribution  for  the  steady  part  is  close  to  that  for  USSAERO 
and  is  thus  acceptable  for  use. 

The  chordvi<=o  AC  distributions  assumed  for  load  trans- 

P 

formation  are  compared  with  the  aerodynamic  distribution  on 
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Table  IV.  Program  Set-Ups 
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Figure  A.  The  constant  distribution  assumed  for  the  equiva¬ 
lent  loading  in  the  first  throe  cases  caused  the  load  to  be 
shifted  aft  as  shown  for  case  3  on  this  figure.  The  linear 
distribution  assumed  for  the  fourth  case  performs  signifi¬ 
cantly  better  in  this  respect.  Similar  changes  in  the  span- 

wise  AC  distributions  are  noted  in  Fiqure  S.  The  linear  dis- 
p 

tribution  for  the  spanwise  direction  transferred  higher  loads 
to  the  coarsely  meshed  Lip  reqion  due  to  extrapolation. 

This  redistribution  of  pressure  over  Lhe  wing  causes 
errors  in  the  transformed  coefficients  of  lift,  pitching 
moment  and  rolling  moment.  To  illustrate  this,  these  co¬ 
efficients  were  calculated  'for  the  four  cases  for  both  aero¬ 
dynamic  and  transformed  forces  and  are  shown  in  Table  V. 
Included  is  the  percentage  difference  induced  by  the  load 
transformation.  The  linear  pressure  distribution  induced 
the  smallest  errors  for  the  load  transformation. 

The  final  deflected  shape  of  the  tip  chords  due  to 
the  flexible  wing  airloads  are  plotted  for  all  the  cases 
on  Figure  6.  All  cases  showed  a  similar  chord  shape,  but 
the  total  displacements  due  to  spanwise  bending  is  differ¬ 
ent.  This  is  caused  by  the  combination  of  aerodynamic  mesh 
size  and  load  transformation  method,  which  created  varied 
wing  center  of  pressure  locations.  It  can  be  seen  that,  for 
this  model,  the  elasticity  efforts  are  related  to  twist  and 
camber  changes. 
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Table  \ 


Comparison  of  Aerodynamic  Coefficient 


1 

CASE** 

1 

CASE 

2 

CASE 

3 

i 

CASE 

4 

COEFFTC 1  ENTS 
DUE  TO 

1  FLEXIBLE 
AIRLOADS 

CL 

. 1324G 

.  13625 

CM 

1 

-.  10995j 

m 

-  .  10596 

L  Cr  . 

n 

. 13508 

.  13932 

■■■■■ 

i  •(  miff  i  i  •  i  ents 

DUE  TO 
TRANSFORMED 
LOADS 

1 

(  ’ 

L 

1 

i 

1  .01  !  JO 

' 

.  i  I  L'i; 

CM 

■■ 

-.09457 

-  .11191 

CR 

.11845 

. 11489 

.  14161 

i 

1 

JM  MHHm 

PERCENTAGE***! 
ERROR  IN 
COEFFICIENTS 
DUE  TO 

LOAD 

TRANSFORMATION 

ERROR  IN 

CL 

-15.74 

-16.94 

1  +  5.54 

ERROR  IN 

CM 

-  7.28 

-  8.60 

j 

1 

j  +  5.32 

ERROR  IN 

CR 

mm 

♦All  the  moments  are  calculated  about  the  origin. 

♦♦Transformed  loads  were  not  available  for  the  sequence 
used  for  this  case. 

♦♦♦For  cases  2,  3  and  4  tho  percentage  error  in  the  co¬ 
efficients  was  calculated  as  follows 

ERROR  =  TRANSFORMED  -  AERODYNAMIC 

TRANSFORMED  X  iU0/' 
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Vertical  Displacement,  Z,  feet 


O - o  Case  1 

A* - ACase  2 

y- - ycase  3 

q — - Q  Case  4 

Y/C= .25 


Percent  of  Chord,  X/C,  Dimensionless 


Vertical  Wing  Tip  Deflections. 


The  solution  convergence  ior  the  four  cases  is  illus¬ 
trated  in  Table  VI.  The  new  DMAP  sequence  used  the  norm 
of  the  incremental  displacement  vector  as  the  convergence 
check,  while  the  change  in  lift  coefficient  was  used  for 
the  interfaced  NASTRAN-USSAERO  sequence.  The  iterative 
process  converged  slightly  faster  for  the  finer  aerodynamic 
mesh  used  for  cases  three  and  four.  The  bilinear  pressure 
distribution  of  case  three.  Each  loop  pass  in  the  iterative 
process  used  in  I  he  new  UMA!’  sequence  required  upprox  i  mu  Lely 
twenty  seconds  of  central  processor  time  and  forty  seconds 
of  input/output  time  for  this  particular  model. 


Table  VI.  Convergence  of  the  Four  Cases, 


LOOP 

CASE  1 
°L 

CASE  2* 
NORM*  * 

CASE  3* 
NORM** 

CASE  4* 
NORM*  * 

1 

.  14192 

18.2335 

16.428 

20.069 

2 

. 13853 

1.627 

.74882 

.92205 

3 

.  13668 

.0788 

.  10017 

.09327 

4 

. 13799 

.0289 

.01048 

.0092 

5 

. 13673 

.0003264 

.00118 

.000837 

6 

.13802 

.0004375 

.000128 

.000094 

7 

. 13762 

.000031 

8 

. 13764 

MAXIMUM 

CHANGE  IN 

I  Z  DISPLACEMENT 
|  FOR  LAST  PASS 
' 

_ 

7 . 81 ( 10-6) 

_ 

-3. 184( 10"6) 

-2.347(10_6) 

♦ELOOP  =  0.0002  for  cases  2,  1  and  4. 

**The  norm  shown  for  each  loop  are  the  incremental  values. 
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VI  Discussion  and  Recommendations 


The  new  DMAP  sequence  has  shown  that  with  some  refine¬ 
ments  it  can  be  a  very  useful  method  of  calculating  flexi¬ 
ble  aircraft  airloads.  The  sequence  could  be  modified  to 
allow  restart  capability  for  large  jobs  and  also  modified 
to  allow  solutions  for  multiple  angles-of-attack.  The 
methods  of  transforming  the  airloads  could  be  improved  for 
quadralateral  plate  elements  by  using  a  bilinear  AC^  dis¬ 
tribution  over  the  element.  This  would  allow  for  the  calcu¬ 
lation  of  consistent  moments  in  the  load  vector  and  would 
better  represent  the  actual  airloads. 

Having  a  module  in  the  program  to  set  up the  required 
matrices  will  make  the  sequence  easier  to  use,  the  data 
blocks  will  not  have  to  be  created  separately  and  then  in¬ 
put  with  the  bulk  data .  This  module  could  be  set  up  to  use 
any  of  the  load  transformation  methods  depending  on  the 
requirements  and  the  model  used.  The  module  could  also  be 
used  to  set  the  initial  angle-of-attack,  and  would  thus 
eliminate  the  need  for  a  real  eigenvalue  analysis. 

The  surface  spline  interpolation  used  in  these  trans¬ 
formations  needs  to  be  looked  at  to  set  how  the  aerodyna¬ 
mic  mesh  can  be  specified  to  make  the  interpolation  of  ACp 
distribution  more  accurate.  The  linear  attachment  flexibi¬ 
lity  parameter  DZ  should  be  checked  also,  relaxation  of 
this  parameter  should  smooth  the  interpolated  surface. 
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The  new  sequence  currently  uses  the  doublet- lattice 
aerodynamic  method.  Although  this  method  can  be  used  for 
transonic  speeds,  the  results  produced  are  not  accurate. 
The  other  aerodynamic  methods  and  the  interference  bodies 
need  to  be  looked  at,  to  see  if  they  can  be  used  in  this 
sequence.  This  would  allow  for  subsonic  modeling  of  the 
whole  airplace  and  the  use  of  wing  models  for  supersonic 
speeds . 
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STATMENT  3  ETJEN3E  R<>R  A  FL  'X.' [ILL  MING  MOO  EL 


LEVE.  2.P  »AST3AN  OiAP  COMPILER  -  SOURCE  LISTING 


OPTIONS  IN  EFF 

1  SEGH 

2 
2 


crt  3 
t 

j  *«*(■*<******. 

t 


E9R=2  NOLI  ST  NOOtC<  NlH£F  NODSOAR 


2 

t 

TMIS  SECTIO  COMMON  TO  RIGID  FORMATS  1  AND  lii 

2 

1 

2 

GP1 

GE3M1  il,jFOM2>  /GPL,ECEyiN,GP0T,CSTM,TGP0r,3U/V, 
Mor.prr  s 

3 

SAVE 

L'JSE"  .NOGPDT  3 

4 

CONO 

ERROR  1 ,NO  SPO  T  1 

0 

GP2 

SEOM7, EOEXIN/ECT  1 

6 

PARA ML 

POO  3/ /C,  N,  PR  ES/C  ,  .V,  /C ,  H ,  /C,  N ,  /  V  ,  N,  JUMOPL  DT  i 

7 

GP3 

GEO  M3,  EOEXIN  ,<"ECM2/,GPTT  /  V,  N , NOGRAV  3 

8 

TA1 

EOT,EPT,OGPOT,SIL,f-PTT  ,C3T«/dST  ,&EI,GPt'CT,/V,r: 
NOSTMP/C, N, 1 /V,K,NOGENL/V,N, GENEL  3 

9 

SAVE 

NOGFNL  ,NOSI»‘  F.GENEl  J, 

in 

CONO 

FRROpl.NOSIMP  $ 

n 

PAKAM 

//C,M, AOD/V, N,NOKGGX/C , N , 1/C ,N , 0  1 

12 

PARAM 

//C  ,M, ADO/V, N,NOMGG  /C,N,i/C,N,0  S 

13 

c  MG 

EST fCGTM,MOT ,OIT , Gr  OK? , /^cL M ,KO ICT, MFLN , HO I C  T  , 
N,  *1 OM".  G  /  C  ,  N ,  /C,N,/r  ,N,/C,Y,COL'PMAE3/C,V,CP*AR/ 
OPT  DAO  1/C  ,  Y  ,  C°0L'6D2/r;f  Y,  CPTRIA1/0,Y  ,CPr-;IA?/C, 
CPOrPLT/C, Y , CPTFPLT /C,Y ,CPTR9SC  J 

14 

SAVE 

MOXF-'X  jNOMSG  E 

1  5 

CONO 

JHPKGGX.NOKGGX  S 

16 

CMA 

GPECT  ,KQICT,  k.*LK,/KGGX,GP3T  s 

17 

LAOEL 

JHPKftOX  $ 

39 


10 

CONO 

ERROR  ijUOMGf.  $ 

19 

EMA 

GPECT,MOICrt  MELM/MGG,/C,N,-l/C,y,HTHAi5:1.0  * 

20 

CONO 

LC°  WG>  GROP.NT  S 

21 

CPWG 

BGPQT,  CSTM,LQEXIN,FGG/OCPWG/V, Y , GRPPNT  =  - 1 'C ,  Y, Mr  MASS  f 

22 

OFP 

OG°WC, ,,,,//y,N,CARONO  S 

2  3 

I.AOCL 

LCPHG  3 

2*i 

t  OUIY 

KGCXjKGG/NQGFNL  $ 

25 

CONO 

L9LU  f  NOGENL  S 

26 

SHA3 

GEI ,/KGGY/V, N,LUSET/V,N, NOGENL /C,N, -1  i 

27 

A  00 

fcCGX.KGC.Y/KGr,  i 

20 

LAPEL 

L3L11  $ 

29 

P  AR  AM 

//C,M,RPY/V, N,NSKIP/C,N,0/C,N,U  3 

30 

PATAM 

//C  ,H,  SU8/V »  N,  Of.  SI  NT /C  >  N  ,8/C,N,  1  ?• 

31 

PUFGE 

BB0BB,OlJE,O2Jt/OESINr  » 

32 

PAR  AM 

//C  |NfMPY/V|NiN0CAf'5ER/CjN>l/C»YjICAMB  =  -i  5 

33 

PURGE 

CAMSER/NOCAMBER  ? 

34 

PARAM 

//C  *N, AOO /V, N, TRANS/C fN, l/C,N,0  3 

3  5 

JUMP 

SETUP  $ 

36 

LABEL 

SETUP  3 

37 

CORO 

LGL41, TRA  NS  S 

3  6  JS*  *■*!•  •*•*****«  •'•■»#*#*•**  »>****  •♦**♦*«#■#*•*•♦**«******•  »#»*♦*  •  ♦  »  ********* #*4* 
38  i 

38  i  THIS  SSCTI3  COMMON  TO  RIGn  FORMAT.'  11 

38  ; 

38  SP't  CA3ECC  ,6F0M%  ,tfU  XlM,GPi)T  ,BCPnT,CSTf'/R-J,  ,U>:r  ,ASU /V,  N,I.US£r/ 

v,K»«!,cn/tf,  u,nrcF?/v,u,  s  ingle/ v,n,  OMjr/Y,N,RSACi  /v ,n,hskip/v, 


40 


STATMENT  SFOJENOE  t'OR  A  FLEXIBLE  WING  MODEL  i 

i 

! 

LEVE.  2.4  SI  ASTRAN  OiAP  COMPILER  -  SOURCE  LISTING  i 


V,REP£AT/V,N  ,NOS ET / V , N, NOL/ V,N , NOA/C , Y ,  SUB  ID  1 


39 

SAVE 

KPC FI, MPCF2, SINGLE, OMIT, R EACT ,NSKIP, REPEAT, NOSE T, NOL, NO A  '$ 

4  0 

GPSP 

GPL,GPST,U3ET,SIL/OGP£T/ V,N,NOGPST  i 

41 

SAVE 

NOGPST  % 

4  2 

CONO 

Ll3L4 ,  N  OGPST  3 

43 

OFP 

OGPST,  ,  , , ,// V,N, CAFONO  i 

44 

L  A0E_ 

L0L4  $ 

4  5 

EOUI V 

tCGG  ,  KNN/MPGF  1/MGG,NNN/MPCF1  S 

46 

CONO 

L9L2,MPCFl  $ 

47 

MCtl 

USET,RG/GM  $ 

48 

MCE2 

U3ET»GM, KGG,  MGG,  ,  /  KliM,  MNM  ,  ,  $ 

49 

LABE. 

L9L  2  S 

5  J 

cOUIV 

KNN,KFF/SINGLE/NNN,«FF/S INGLE  i 

31 

CONO 

L3L3, SINGLE  ! 

52 

SCE1 

US  ST. *  NN , MNM  ,,/KFF,KFS, ,MFF , ,  % 

53 

LABEL 

L3L3  3 

54 

EOUIV 

KFF ,KA A/OMIT /  MFF,MAA/OMIT  l 

55 

PURGE 

G3/0'*IT  3 

56 

CONO 

L1L5,0MIT  $ 

5  7 

PARAM 

//fi  ,*>,  PRSC/ V  ,N, PREC  J 

58 

VEC 

U3ET/V/C,M,F/C,N,0/C,N, A  I 

59 

PAUX 

»V»/KOO,  ,  KOAfKAAB  i 

5  L 

OECOIP 

KOO/t  00,1)00  /  C,N,  1/C  ,  N,0/  V  »N  ,  N1NI1/V,  N  »DET  /)  ,  V  ,  NOE  f  /V  ,N,SING 

51 

SAVE 

MIND, OET , NOE  T ,51 NG  i 

62 

FBS 

LOO  ,UOO,KOA/r,0/C  ,H»  1/C  ,  N  ,  -1/ V,  M  ,  PREC/V,  U,  PREC  5 
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LEVE.  2.1  ‘I  AST  RAN  DMA  n  COMPILER  -  SOURCE  LISTING 


53 

HPYA3 

6  4 

SMP2 

65 

LABEL 

5  6 

COLD 

67 

RBEG1 

60 

KUMI.2 

59 

RBHG3 

70 

RRE'G'+ 

71 

LABEL 

72 

DPO 

73 

SAVE 

74 

CO  NO 

75 

cnuiv 

76 

REAO 

77 

SAVE 

7a 

OFP 

79 

COLO 

ac 

N  T  RXI H 

ai 

SAVE 

82 

touiv 

83 

G  KAO 

04 

GKAM 

K3A,G0,KAAB/KAA/C,L,1/C,N,1/C,N,1/V,N,PREC  S 
USET,GO,HFF/ MAA  S 
LBL5  $ 

L3L6, REACT  S 

USET.KAA ,MAA  /KLL , KLK , KRK , MLL ,M L R,MRR  J 

<LL  /I  LL/  i 

LLL  jKLR.KRP/OH  i 

DM, HLL  »MLR, MFR/MR  £ 

L0L6  S 

DYNAMICS, G&L  ,SIl , USE! /G^LO, SILO, USFTD, TFP DDL  , ,  ,  ,  ,  , E ED, EQD YN /V , 

M. LOSET/V  ,N,  LMSETf  VV,L,  MOTFL/V,  N.NOOl  T/V,M,  NOFCDl  /V  ,  N,  KlOPRL  /V  , 

N,  NOW  FT  /  V ,  N  ,NC1  .\L / V , U,  HOLE O/C  ,  N , It  3/V,N,N0'JE  » 

LUSETO,NOUE, NOEf 0  1 

ERRORS  ,’IOEEO  t 

GO, GDD/NOUE/firt, GMO/NOUE  i 

KAA,NAA,MR,OH,CE.O,USET,  CASECC/LAMA,  PHIA,  *1  ,  OSIGS/C,  L,HOOES/V,N, 
NEIGV  3 

NSIGV  S 

DEIGSjLAHA,  ,  ,  ,//V,N,CARCMO  i 
ERRORS, NEIGV  $ 

CAS  E(CC  ,MA  TPO  OL  ,  C  Ol'YfJ ,  ,  rFP0UL/K2PP,|1  2P»VJ?I’P /V  ,  M  ,  L  US  ETD/ V  ,  N, 
NO<SPP/V,N,NOH2PP/V,N,NOS2PP  S 

N3K2PP,  NO  M. 2  PF, 110  02H*  S 

M2PP,H200/N0  A/C2PP,Dcl>D/MOA/K2PP,,<2fJO/tnA  i 

USTTO.GM, GO,  , ,  ,,K?PP,M2PP,B2PP/,  , ,  S  DT  ,  <  ?  Ci",  k  ?P0  ,  ^200 /C  ,  M, 

CM  =  LiV/S,M,  n  Iff. /C,N,  *100/1.  /.;  ,  N,  0 . 0/C  ,  N,0  .i  t  '  ,  L  ,0 .  O/V  ,N,  NO*  ?X>  P/  V , 
N,*»nK?PP/V,N  jMCftCPF  /  V,L,rT>CFl/V,N,  SINS*?  /  V  ,  N  ,  UM  IT  /  V,H  ,  NOUfc/C  , 

c,n,-i/c,u, -i  i 

U3ET0,  PH  IA,,  LAMA,  Dll  , C2 0 0, M 2 00, K2PD, CAS tCC/MHH,BHM,  KHH,  AAA  A  A/  V, 
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5T  ATMcNT  S  E0JEN3E  r0K  A  FLEXIBLE  WING  MODEL 


LEVEL  2.3  W4STRAN  OMA  v  COMPILER  -  SOURCE  LISTING 


N,  NO'»r /C,  V,L  MOCES  =  O/C,Y,LFREO=0./C,  Y  ,HFR£a=0./‘/,  N,  NOM2PP/V,N, 
NOB?PP/V,N,NOK2FP/V,N,NONCtJP/V,  N,FMODE/C  ,  i  ,KOAMP  i 


8  li 

SAVE 

NONCUP  »FMOOE 

86 

A  00 

A4  AAA,  CAMBER  /PHI OM/C , Y , A  0 AP= <1 . G , 6.  0)/G,r,nETAr  (1,0,C.'))  $ 

87 

APD 

F9T  ,6CLDYM, EOT, BGPOT, SILO  ,  US  ETO  ,  CSTrt  ,  GPLO  XE04  El- 0 ,  E  CT  A,fl6PA  ,S  JLA, 
USETA , SPLINr ,AEfiO, ACPT , FL1GT  ,CSTMA,GFLA ,«  i LGA/V, H ,N Y/V , N,NJ/V, 
N, Ln^ETA/V,N  ,90V  1 

38 

SAVE 

NK, Ml, t  USETA ,90  V  l 

89 

PARAM 

//C,N,MPY/V, N,PFILF/C,N, U/C,H,1  S 

9  C 

COf.O 

S<P°1  T  ,JIIMPPLOT  % 

91 

PARAM 

//C  , N,M°Y /  V,  MjPLT KLG/C  ,  N., 0/C ,N ,  1  S 

92 

PLTSET 

°C0  8,  EQA“R3 ,  EOT  A/PLTSET  A  ,  PLTPAR.A,GPSCTS  A  ,  EI.S  ;TS  1  A  V,  N,NSIL1/9,N 
J(JMPrLOT  $ 

93 

SAVE 

NSTL1  ,  JUMPPLOT  $ 

94 

PRTMSG 

PLTS"TA  //  i 

95 

CORD 

S<PP|.T,JUMPPLOT  $ 

96 

PLOT 

PLT  PARA, GOSETSA , EL  SETS A , CAS  ECC , OG°A , ESA  HPO ,  , , , , / PL CTX  2 /V ,N , 
N5IL1/V,\I,LUSETA/V,U,JUHPPL0T/V,N,BLTFLGAV,N,PFILE  1 

97 

SAVE 

PFTLF, JUMPPL CT,PLTFLG  3 

98 

PPTMSG 

PL0TX2  //  i 

9  9 

L  APE.. 

S<°PLT  3 

18  0 

COKO 

t RR0R2 , NOSED  S 

19  1 

GI 

SPLTliE,USir  ,CSIMA,0GRA,S1L  ,  ,GM,60/Gf  KA/V,  M,N</V,  N,L'JSET  5 

10  2 

P  ARAM 

/f C,M.  Am/V,N,DtSJRY/C,H,Q/C,N,l/  l 

13  3 

AMG 

AERO,ACpF/AJJL,SKJ,  01 JK  ,  02  JK/V  ,  II ,  ‘l*/V,  H,  N  1  /  V  ,  M,  P  E  1ST  PY/  £ 

1>J  h 

SAVE 

OESTRY  S 

10  5 

PARAM 

7/C,f»,SU0/V,  N,IHEO/C,M,0/C,N,1  S 

186 

3 
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LEVE.  2.0  NAS’RAM  DWAP  COMPILf.f  -  '.OUWCt  LIs.'ING 

10  6  S 

10  6  I 

106  'I  THE  KOQULE  WILL  REPLACE  THE  FOLLOWING  sTAMENTS 

10  6  S 

106  I  T'i  C  PH  AP  INST  r-‘UC  T 1 C  N  f'jk  THE  NOODLE  WILL  JE 

106  3 

106  T  Gt  AS  ArPT,T03eL,SCTA,GFLA ,  TOP  A /S  AS  ,  GT  AK ,  7GLS  ,  V  CMS ,  VTM  A,  VT  RA  , 

106  S  VTM,VTR,AIOMIT/5 

10  6  : 

106  EQUIY  SEE  ,  3A S/IHEO  5 

10  7  T RNS3  GTAXT/GTAK  2 

io  a  s 

108  i  THIS  ENOS  THE  SCCTION  TO  PE  REPLACED  ftf  THE  MODULE 

10  8  t 

136  ».**¥»*»*****♦•*  *►♦***#•**»*****»♦».»  .•*  ♦****.«♦*  .»*♦*».  **.♦»»«»»¥.*«  *****♦»*. 

10  6  S 

108  EQUI7  OlJXjCPPH/OESINT  i 

109  A  00  SAS,CP°M/TMA /C,Y , ALPNA= < 1 .8 , 0. 0 ) /C, Y ,PET A= ( -1 .0,0.0)  i 

110  MPYAD  TH»,5KJ,/TrAA/C,N,l/C,.N,l/C,N,0/C,N,l  i 

111  MPYAD  TfTTT.SASEL, /TT5AS/C.,n,0/C,N,1/C.W,0/0,S,1  b 

112  MPYAO  TTSAS, GTAK,7TSTAP/C,H, 0/0,9, 1/C, N, O/C, M,  1  3 

113  ADD  T3TAF,B9nDO/7STA/C,V,ac/C,Y,0tTA=<6, 0,6.1)  i 

114  MPYAD  AI0MT,TTAA,//VTNA/C,N,1/C,N,1/C,M,0/C,N,1  > 

lib  MPYAD  TTAA,7TKA,/T06AA/C,N,0/0,N,l/i;,N,0/C,N,i  i 
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ST  A THE  NT  SrOJCNIF!  FOR  A  F|.  CXIBLC  W]!.G  MO1)  EL 


LEVEL  2 •  r  NAS'RAN  DMA  F  COMPILER  -  SOURCE  LISTING 


11 0 

MPYAD 

TEA  A, VTMA./T VEMA/C , N , 0/C , N, 1/C  ,  N  ,0 /C  , N ,  1  » 

117 

PAP  AM 

//C ,N, A09/V, N,XQHHl/C,N, 1/C,N,0  $ 

116 

AMP 

A JJL,SAS,D1 JK,02JK,GTKA, =HIDH,01 JE, 02  J  E  ,  J  SE  T  0  ,  A  E  RO/ OHHL I , 
QH\JH/V,N,NOUE/V,N,  XQMIIL/V,  Y,GUSTAERO  =  H  e 

OKHLl , 

119 

SAVE 

XIH^IL  $ 

4 

12  1 

MPYAO 

C°P^,QKHLI,/  <1KKXI/C,M,1/C,N,1/C  ,II,»/C,N,  1  1 

12  1 

MPYAD 

TS  T  a\qKKK I ,  /  PC/C  ,N,  0/C,N,  i/C,N,  0/C,H,l  I 

12  2 

PAR  A  MR 

//C  ,N,M°Y/V,N,RCAF.IA/V,Y,RC/V,Y,WAREA  $ 

12  3 

PARAMR 

//C  ,N, MPY/V, N,OAREA/V,Y, MAREA/V,Y,a  $ 

12  »» 

PARAMR 

//C  ,N,  MPY/V,  N,ORCAR.£A/V,N,KCAREA/V,Y,Q  5 

12  5 

PARAM 

//C  ,N,  MPY/V,  N,TRANS/C,ll,l/C,N,-1  S 

12  6 

COMO 

L9L10, TRANS  J 

12  7 

LAGEL 

L3LMJ  $ 

12  8 

»*¥***  .»•»«****  »♦*»*»»-*•-■.***♦«•>♦»»*♦«*»*»*«  |  ,  *  »  »  .  «.  •  M,«<)  »** 

12  6 

s 

126 

$ 

THIS  SECTXO  COMMON  TO  RIGID  FORMA'S  1 

12  6 

< 

.*■ 

12  0 

SPA 

CASFCCjGEOM*  ,  EOEXIN,  <"-PDT  ,  99  POT  ,  C'JTN/RG,  YS, US  cT  ,  A  S  ET  /V, :)  ,LUS  FT / 
V,N,h|PCFl/V,  S,KPCFr/V,N,  SINGLE/ V  ,  N,  OMTT  /  V  ,  N  ,  AE  AC  l  /V  ,N,NS*IP/V, 
N,  RET  EAT/V,N  ,NOS  ET  /V,  K,  MOL/  V ,  N  ,  NOA/C  ,  Y ,  StmO  i 

129 

SAVE 

M°C  FI , MPCF2 , SINGLE, OMIT, RLACT, NS<T P, RE  =  t U ,NGSrr, MOL, MO A 

i 

13  vi 

CONO 

ERR0R3  ,NOL  3 

131 

PARAK 

//C ,N, AHO /V , N, NGSR/ V, N, S INGLE/ V,N, RE ACT  { 

132 

PURGE 

KRP  jK.LRjOR.OM/REACT/GV/mpcFI/GO  ,  KCO  ,  LC;3,P0,U00  7,  *  UP  V/O'II  T /P  5, 
KPS, IC3S/SING  Lt'/OG/NOSR  1 

13  3 

COMO 

LDL1A,GL‘IEL  S 

13  A 

GPSP 

GPL ,GPST , USE T,SIL/GIFST/V,M,NOGPST  5 

Arj 


sr  a  tug  n i  seouencg  r or  a  flexible  wing  model 


LEVEL  2.0  NASTRAN  DMA  P  COMPILER  -  SOURCE  LISTING 


135 

SAVE 

N3GPST  S 

136 

CONO 

LOLlf> NOGPST  1 

137 

OFP 

OGPST,  ,,,,//  3 

13  8 

LABEL 

LBL14  S 

13  9 

EOUIV 

KGG  » KNN/MPCF 1  i 

14  0 

CONO 

LtlL 15  r  MPC F2  3 

141 

MCEi 

USET,RG/GM  $ 

14  2 

MCE2 

USGT.OH.KGG,  ,,/KNN,,,  $ 

14  3 

LABEL 

LBLlb  3 

144 

EOUIV 

KNN  »  KF  F/S IMG  LG  3 

14  5 

COhO 

L1L 16 .SINGLE  3 

14  b 

SCE1 

USET  ,KNN.  ,,  /  Kf  F  ,KFS,  ;SS  ,  ,  , 

147 

LABcL 

LOL 16  3 

14  8 

EOUIV 

kff.kaa/omit  « 

14  9 

CONO 

LPL17.0MIT  5 

15  3 

SMP1 

USET,KFF» i,/GO,KAA .koo.loo 

151 

LAGtL 

LOL  IT  3 

15  2 

EOUIV 

KAA  .KLL/REACT  1 

15  3 

CONO 

L9L16, REACT  3 

15  4 

RBMG1 

USET.KAA, /KLL.KLR.KPR,,  ,  4 

15  5 

LABEL 

L0L16  i 

15  6 

R8MG2 

KLL/ILL  t 

15  7 

CONO 

LBL If. REACT  3 

15  8 

RBMG3 

LLL,KLR,KR*/OM  S 

15  9 

LABEL 

LBL 19  S 

46 
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STATMLN1  SETJENCE  '  OK  A  FL  EX19LE  WING  MO!) EL 


LEVEL  2.0  NASTRAN  DMA P  CO) FILER  -  SOURCE  LISTING 


SETUP, i  < 

PC,/PGI/0£SINT  £ 

ULVT,t)OOVI/I  HEQ  5 
PSI  ,QRI/IWEQ  ? 

ITERATIVE  LOOP  FOR  FLEXIBLE  WING  CALCULATIONS 

LOOPTOP  3 
LOOPTOP  3 
FINIS, ELOOP  ? 

USET,GM,Y.S,KFS,GO,DM,PG/ftR,PO>PS>PL  S 

LLL  ,KLL,°L,L  CO,  KOO  ,  PC /PHI OH  ,  UOO  V  ,RUL  V , R'JOV/ V  ,  N ,  OM IT/ 
V,Y,IRES  =  -1/ V,N,NSKIP/V,N,EPSI  ‘i 

EPSI  $ 

LBL20 , IRES  $ 

GPL ,  USET , SIL  ,RUL V//C , W, L  S 
G°L ,USET ,SIL ,RU0V//C,N,0  i 
LI3LZ0  3 

PHTOP,ULVI/ULVT  $ 

PS,PSr/PSTT  s 
UO0V,  JOOVI/WOOV1  2 
L0L3O, REACT  J 
aRjftRI/'SRT  < 

LOL30  $ 


ISO 

LEVEL  ; 

KEPT 

161 

COPY 

16  2 

PURGE 

16  3 

PURGE 

164 

t 

154 

<i,  *****  i 

15  4 

S 

164 

s 

16  4 

$ 

164 

JUNP 

15  8 

LAPEL 

156 

COND 

157 

SSG2 

168 

SSG3 

159 

SAVE 

17  U 

COND 

171 

UATG“R 

172 

HATGPR 

17  3 

LABEL 

174 

A  00 

175 

A  00 

17  6 

AOO 

177 

COHO 

17  8 

AOO 

17  9 

LAlUu 

STATMENT  SEOJENCE  TOR  \  Ft.  EX1BLL  Ml  NO  MuOEL 


LEVEL  2.0  94S'RAN  DMA"5  COMPILER  -  SOURCE  LISTING 


tao 

PARAM 

//C  ,N,  AOO/V,  N,X0HHL/C,U,1/C,N,G  j> 

101 

AMP 

4 JJL,  SAS,  91 JK.02 JK,CTKA , PHI  OH, 0 1 J£, P2 JE , USET 0 , A ERC/ OHHL , QKHL, 
OHJL/V ,N,NOU£/V,N, XOHHL/V, Y, GUST AERO  S 

18  2 

SAVE 

XQHHL  $ 

18  3 

MPYAO 

CPPM,QKML,/aKKK/C,N,i/C,N,i/C,N,0/C,N;l  i 

18  0 

M  PYA9 

TSTA,ClKK<<,/PG/C,U,0/C,M,  1 /C  ,  N,  0 /C,N  ,  1  $ 

18  6 

A  00 

Pfr, PGI/PGT  $ 

186 

A  00 

OKKK, RKKKT/QKKKT  3 

13  7 

MPYA9 

PM I D-l i  DHIDII,  /7N0R V/C , N ,  1  /C ,  N ,  1 /U  ,H,  0/C  ,N ,  l  i 

18  6 

P ARAML 

TNQRV//C,N,OMI/C , N, 1/C, K ,1/ V , H , S«MN  % 

139 

PAKA.MR 

//C ,N,SftRT/V ,N,VNSH/V,N,SOVN  $ 

19  h 

PRTPAWt 

// C VMSR  3 

191 

PARAMR 

//C  ,M,  LT/V,M  »CHL  f  OF/ V, II »  V  TSR/C  »  Y  ,CL  OOP//  / /  V  >  N ,  F L  A G  i 

19  2 

CONO 

L0L50,  FLAG  i 

19  3 

PARAM 

//C  , ft,  AOO/V,  N,IWEa/C,H,0/C,N,l  3 

19  4 

PURGE 

ULVT, PSI .UOOVI/IHEQ  2 

195 

CONO 

LBL3i,HtACr  S 

19  6 

PURGE 

RRI/1MEQ  3 

197 

COPY 

ORT/QRI/OESI  NT  t 

19  8 

LAPEL 

L8L31  S 

199 

COPY 

QKKXT/OKXKI/ OESINT  1 

2U0 

COPY 

l/L V T/ULVI /  Ot SINT  3 

201 

COPY 

OGT/PGI/OESI NT  t 

20  2 

COPY 

PST  T/PSI / OES  I NT  3 

21  3 

COPY 

U90 VT/UOO VI/ OESI NT  * 

2U4 

REPT 

LOO  PTOP,  10  S 

48 


STATMENT  SEQUENCE  r0R  A  FLSXIGLt  WING  MODEL 

L  EVE.  2.0  9ASTRAN  OMA®  COMPILE?  -  SOURCE  LISTING 

2U5  LAOEL  L3L50  5 

20  6  *  4-  **  >  *  r-  *9  **  *  ******  *■**  »  *#  *  **  V  **  *  «*  ***  *  ***■  *  *  *  *  #  *  »  *  *  #  *  *  *  *  *  *  v  *#  *  *  »*  4*  *V  *  * 

20  6  $ 

236  !  TS1IS  SECTION  SIMILAR  TO  RIGID  FORMAT  1 

206  J  FOR  THE  RECOVERY  CF  TOTAL  OISPL AC EME.'I T S  A  NO  STRESSES 

20  6  l 

236  MATPRN  OKK XT ,  Ul. V T ,  ,  ,//  % 

2D7  PARAM  //C,N,A00/V,N,NSKIP/C,N,0/C,N,1  2 

208  JUMP  RECOVERY  i 

209  LAOEL  RECOVERY  X 

2U  SOR1  USET,PGT,ULVT,U0OVT,  YS,GO,GM,PSTT,Kr  S,ICS$.QRT/UGV,PGG,aG/ 

V,N,NSKIP/C, N, STATICS  3 

211  PARAM  //C,N,AOO/V,N,NSKIf»/C,N,l/C,N,l  S 

212  REFT  RECOVERY, 1  3 

213  GPFO*  CASECC,UGV,KFLM,  KOICT ,  EC T  ,EOEX  I  N  ,G?£CT  ,P  66  ,QG/GNRGY  1 , OGPFII1  / 

C,N, STATICS  5 

214  OFP  OMRGVlfOGPF91,,,,//  S, 

215  SOR2  CASFCC.CSTM,  HPT, HIT, CttEX IN, SIL ,GPTT, LPT, ?GP0T, ,QG , UGV, SST , 

X  Y  C  00,  PGG/3PGl,0(TGl,0UGVi,0ESl,0EPl  ,  FUG  V 1  fZ  ,  N»  S I  AT  I CS/ V  ,  N  , 
NOSO'T  2=”1  '! 

216  SAVE  NOSORT 2  % 

21 7  COMO  L9L  2! , NOSORT  2  i 

218  S0R3  JOG VJ , OPG1, 0 QG1 , OE FI ,OES 1 ,/OUG V2 ,OPG2 ,#0 J2 ,0 EF2 , 0 FS 2 ,  S 

219  OFP  DUG V2,  0°C2, 0 RG2,OLF2, OES2,//V,M, CAR  PI  10  j, 

22 C  SAVE  CAROHO  f 

221  X  YT  RAN  XlfC  n»,0DG?,  0  0G2  ,  QUGV2 , 0  C.S  2,  OEF  2 /XY=  L  TT /C  ,  M  ,  T  Rf  U /C  ,H  ,PSP  T  /  V  ,K, 

PFILE7V,N,CA RUNO  5 

222  SAVE  PPILS.CARDSO  J 


Si'  ATM  EMI  s  E'l  J  END  E  r0R  A  FLEXIBLE  WI/sG  M09EL 


LEVE.  2.0  NAS" RAN  'VIA  P  COMPILER  -  SOURCE  LISTING 


22  3 

XYPLDT 

XYPLTT  //  S 

224 

JUMP 

0°L  OT  $ 

22  5 

LABEL 

L9L21  3 

226 

OFP 

0«5 VI  ,OPG1,OQG1,OGF1,OE51,//V,U,CARONO  i 

227 

SAVE 

CARONO  $ 

22  0 

COMO 

P2,  JVM  PPL  or  1 

22  9 

LABEL 

t^LOT  J 

23  b 

PLOT 

DlTPAR,GP<5ETS,ElSETSfCASECC,BGP0T,taEXIN,SIL  ,  FUGV  1,  ,GPECT,06Si/ 
PLOT** «./’ V  |N(NSIL/VjN>LUStT/V»Ni  JUilPPL  OTT  t  v ,  N ,  P  LT  FL  5 /V  ,N,PFILE  I 

231 

SAVE 

PFILt  3 

232 

PRTMSG 

PLOT*?//  S 

23  3 

LAPEL 

P2  3 

234 

HUM! 

234 

t 

234 

$ 

THIS  SEDTIDN  CALCULATES  THE  AEROOYHAMID  CO  E  F  F I D  t  E  NT  S 

234 

s 

234 

MPYAD 

ATMATaKKKT,/TFACP/C,N»0/C,N,l/C,M,0/CfN,  i  I 

235 

PAk’AML 

TFACP//C, NiDN:/C,N,  l/C,N,l/V,N,  CLAA  3 

236 

PARAMR. 

//C  ,N,  OIV/V,fl,  CEAGP/V,N>CIAA/V,  Y,WAREA  v 

237 

MPYAD 

TVEHA, QKKK r,  /TMACP/C,N,l/C,N,l/C,N,0/C,Nf  l  * 

23  6 

MPYAD 

TVS RA»  OfcKKT , /TRACP/C ,N, 1 /G, N , 1/C  ,N,0/C ,N , 1  S 

239 

PARAML 

riArp//c,N,om/c,N,i/c,N,i/v,N» acma  $ 

24  0 

PARAML 

TRACP//C,N,nm/C,N,l/C,N,l/V,N,AC*4  i 

241 

PAEAMR 

//C,N,OIV/V,  N,CMACP/V,N,  ACM  A/V  ,  IJ ,  PC  AR*A  i 

24  2 

PAR  AMR 

//C  ,N,OIV/V, N,CRACF/V,H, ACRA /V , N ,F C ARE  4  l 

24  3 

MPYAD 

VCLS,P«T,/TFSGP/C,M,1/:,H,1/C,|I,0/C,N.1  ? 

ST  ATM1. NT  SENTENCE  TOR  A  FL  EXIOI.F  WHIG  HOOEL 


LEVEL  2.0  ’4 AS  RAN  OHA  P  t.  )MP1L:R  -  SOURCE  LISTING 


24  4 

PARAML 

T PS  CP/  /c  ,  N,  0  MI/C  ,  ti ,  1/C  ,  N  ,  1/  V ,N  ,  ACL  5  3 

245 

PAP.AMR 

//C,M, OIV/V, N,CLSGP/V,N, ACLS/V ,  N,QAREA  i 

246 

MPYAO 

VGMS,PGT,/FMSGP/C,N,l/C,N,i/C>N,0/C,N,l  i 

247 

I1PYA0 

VTM,FMSGO,/TNSGP/C,N,i/C,N,  1/C ,  N  ,0/C  ,1*,  l  i 

24  0 

A1PYA0 

VTR,FMSSP,/TRSGP/C,N,1/C,N, 1/C , N,0/C ,N,1  { 

24  9 

PAKAMl 

TMS6P//C,N,DMI/C,N,l/C,N,l/V,N,r}ACMS  3 

25  0 

PARA ML 

rRSGrl//C,NtONI/C,N,i/G,U,l/V,N,OACRS  3 

251 

PAR AMR 

//C,N,  OIV/V  ,  N,CMSGP/V,U,OACNS/V,N,ORCMUA  5 

25  2 

PARAMR 

//C ,N, OIV/V, N, CRSGP /V »  N  ,  0 AGRS/ V , M, ORCARL  A  » 

25  3 

PRTP4RM 

//C ,N,0/C,N, CLACP  J 

25  4 

PKTPARfl 

//C ,M, 0/C ,N , CHACP  5 

25  5 

PRTPIRM 

//C ,N,Q/C,N, CRAC  P  S 

?5  6 

PRTPAKM 

//C ,N,0/C,N, CLSGP  i 

25  7 

PRTPARM 

//C,N,0/0,N, CNSGP  J 

2j  3 

PR7PARN 

//C,N, O/O.N, CRSGP  } 

25  9 

JUHP 

FINIS  3 

254 

LABEL 

ERROR?  $ 

251 

PRTPARM 

//O  ,N,  -2/C,N  » FLUTTER  i 

25  2 

LABEL 

FRRORt  3 

25  3 

PR1PARM 

//S,M, -1/C, H, FLUTTER  i 

254 

LABEL 

ERR0R4  3 

255 

PHTPARM 

//C ,M, -4/C, H, FLUTTER  t 

256 

LAC-El 

ERRORS  S 

257 

PRTPARM 

//O ,H, -3/C, N, STATICS  * 

25  0 

L  AOtL 

FINI$  $ 

51 
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52 


APPENDIX  B 


The  structural  model  used  to  test  the  DMAP  sequence 
consisted  of  rods  CROD,  quadralateral  membranes  CQMEM2, 
and  shear  panel  CSHEAR  elements.  The  CROD's  provide  ten¬ 
sion,  compression  and  torsional  stiffness.  The  CQMEM2 
elements  have  a  finite  inplane  stiffness  but  do  not  resist 
bending  and  the  CSHEAR  elements  resist  tangential  forces, 
but  do  not  resist  normal  forces.  The  elements  used  in 
this  model  do  not  handle  rotational  degrees  of  freedom, 
so  these  were  constrained  out  of  the  solution.  The  number 
ed  grid  points  are  shown  in  Figure  7,  illustrating  the 
relative  shape  of  the  model.  The  elements  were  shrunk  by 
50  percent  to  show  their  relavtive  positions  in  Figure  8. 

The  aerodynamic  mesh  size  used  for  the  panel  with 
twenty  boxes  is  shown  in  Figure  9  versus  the  projected 
mesh  of  the  structural  model.  The  sixty-three  box  panel 
is  shown  in  Figure  10.  The  aerodynamic  models  used  the 
reflected  wing  in  the  calculation  of  aerodynamic  forces. 


Depiction  of  Element  Connection. 


56 


Figure  9.  Projected  View  of  Structural  and  20-Box  Aerodynamic  Model. 


< 


Structure 

Aerodynamic 
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SU^P  jUTI N  £  GIAS 

INTEGEF  it  S, GT AK,SJ SEL ,T71TT,  VGl  S,VGMS,VT*A,  VT“, VTB, VT R A ,  AIOMIT 
I  NTE3EP  A  CPT,TDSE  ,  £DT  A ,  0  *LA  ,  3GFA  ,  PUF 1 ,BUF2  ,  BUr3  ,  CORS  7 
I NTLGtF  TRAIL*  H  r' A  *  1C  3D , TYPI N, T YFOUT 

DIMENSION  TPAILt?)  ,Hrj3A(7>,KC3B(?)  ,7(i>  , 17 ( 1 > , NA HE< 2 > 
equivalence  (rT(i’,7(D) 

COMMON  //'■. 

COMMON  7  S  Y  5  TTi/  T'l,c 

COMMON  AP*,  CKV/TY~TN,7YOOUT,IR0W,NRQW,NM 
N*=l 

DATA  N/ME/'..  HGIAS  wp”5D/ 

DATA  AC3T,TOSEL,-CTJ,GPL  A,BGPh/1^  1, 1  J2,  1L3, 1 JA,  US/ 

DATA  SAS,3T  AK,T"n"r7,  SAS-  l,VGLS,  VG^S,  VTMA,VTRA,\/TM,  VTR,  AIDMITX 
♦  2  o  1, 2,2,233,  2  )  *4  ,  2 1  ■? ,  20  &  ,2L?,2i&,0J9, 211,211/ 

L  CORE  =C0RST ( 7, A ) 

0UFl  =  Lro<E-IP'Jc  +  i 
IF  (BUFl.l E.">  3D  TO  ? 

AN  “N**  15  BEFORE  'ME  DATA  BLOCK  NAME  SPECIFIES  THE  FIRST 
AOOFESS  DF  THE  3L'C<  IN  0  FEN  COF  E 

AN  "L"  BSFORF  the  DATA  3LCCK  NAME  SPECIFIES  A  PARAMETER 
WHICH  IS  T  -»E  LEN-'H  D  F  tm*  PATa  BLOCK  IN  OPEN  COKE 

N  ACPT  =  1 

CALL  GOPEN(  ACPT,  7(PUFi)  T 
CALL  READ  (ACPr,7  (•  AC^T)  ,F  ,1,M) 

N I  IS  HE  NU^ER  OF  A£*0  BOXES  AND  POINTS 
N 1=7 ( NAC3T*tI 
CALL  CLOS  p  l  AC°T , i) 

TRAIL  (1)  =  T05CL 
CALL  RDTR  L ( T °8 IL ) 

ND  IS  THE  HUH3EP  OF  STRU«*TU>AL  £ LEHENTS  LOADED  AND  ALSO 
THE  NUMBEF  OF  ’EE°EN3ENT  STRUCTUFAL  POINTS 

ND=TRAIl(2> 

NTOSEL=NACPT 

CALL  GOPEN(TOSEL»'Tf PUci)  ,  0) 

CALL  READ  ( T OSEL ^ CNTOSEL * ,ND , 1 , M) 

CALL  CL05ECT0SEL  *1) 

T KAIL (1)  =  G3L  A 
CALL  RCTRL(TOAIL) 

LGPLA=TRAILC  2) 

L  BGPA  =L*L  GPL  A 

CALL  GGPENC5PLA,7<B«F1),0> 

NGPLA  =  NT3SEL  ♦NO 

CALL  REAO(G°LA,TC‘rGPLA)fL  GPLA,1,H> 

CALL  CLOS  E  (  SPL  A  ,  i) 

N8GP#  =NG3LA*1_G°L* 

CALL  G0PEN(BGP\,-'<DIJF1),  '  ) 

CALL  READ  (3SPA  ,  7  (  N9G®  A)  ,  L  EGPA  ,i,M) 

CALL  CLOS  E  <  BGD  A,  1"  ■ 

NCORE  =  BU" i-l 
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o  n 


CHECf  TO  SEE  Ip  ENO'JGH  0°^  CORE  IS  AVAILABLE 
IF  (NC0RE.LT.(NBGP1+LTGPA  ))  GO  TO  SjC 
I  CONE  M2*5  3  3  IS  TuE  INTERNAL  ADDRESS  OF  CQMEM2  EL 
THE  FOLLOWING  S  S’ATMENTS  LOCATE  THE  RECORD  FOR  T 
ELEMENTS  IN  THE  E‘TA  TA3LE 
CALL  GOPEK(ECTA,^»aun),  ) 

IC0M£M2=r3  3 
N  ECT  ft  =  N9S  FA  «■  LBS3  A 

1  CALL  READ  (ECTA,T/. NECTA)  ,  i  ,  M) 

IF  (I Z ( ME  CIA  1 • EG .1C CM" *2)  GO  TO  2 
CALL  RtAD(ECTfi,7eNECTft)  ,l,l,M) 

GO  TO  i 

2  TRAIL  <1  )s$IL 
CALL  FDTRL  (T  RAIL  ) 

NDOF=TRAIL(?) 

NGP=ND0F/6 

BUF2=PUF1-I:»tic 

BUF3=QUF2-IBM“ 

N CORE  =  PU-  3-1 
L  ECT  A  =  7 

NNC=NECTA  +LECTA 
L  ND=2  *  NO 
NTEX=NNO*LNO 
LTEX  =  NGP 
NTEY  =  NT  EX+LT  EX 

ltey=ngp 
NSASE L*LT  EY^MTEY 
LSASEL  =  ND 

NTTTTT=NSf SEL+LSASEL 
LTTTTT  =  NOCp 

nvgls=uttttt*lttttt 

L  VGLS  =  NOOF 
NVG‘'S  =  NVSLSM.VGLS 
L VGHS=N03F‘HGO 
ICOREN=NVGTS+LVG«S 
IF(NCOF.E.  L7.IC0RFM)  GO  TO  GOl- 
FLAG=-1 C  .  C 
DO  9  J=l,  NODF 
7  CNTEX+J-l)  =FLAG 
7CNTEY  +  J-I)=FLA5 
9  CONTINUE 

00  1C  J=1 » L VGH5 
Z  CNV5MS ♦J-l) =n. 3 
1C  CONTINUE 

DO  7  J=l»  LVGLS 
ZCNVSLS+J-i>=0.1 
7  CONTINUE 

CALL  G0PEN(SA3EL,^(9UF2)  ,  1) 

CALL  GOPENCTTTTT  >Z.(9Ue3)  ,  1) 

TTPIN  =  1 

TYP0UT=1 

IR0W=1 
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T  lil 

X 


MCqA  ( 1)  =S  ASE L 
MC8A(2>=* 

MCSA  (  3 )  =L  SA  3  f L 
M CHA  C  **  >  =3 
MCPA(E)=t 
MC8A  ( 6  )  =. 

MC°A  (7 )  =  . 

MCBn ( 1 ) =TTTTT 
M  C  B8  (  2 )  = ' 

MCt?8(3)  =L  TJ  TTT 
MCB8<4>=2 
MCBP<‘  )=1 
MC33(6>=. 

MCBB(7)=: 

LOOP  3  CALCULATED  TV“  FOLLOWING  DATA  BLOCKS  SAS 
TTTTT  »  VGLS*  A  WO  VC'S 
DO  3  1  =  1,  f  O 

4  CALL  P.EAQ(£CTA,7f  ’FCTA)  ,"•  ,T,P) 

IFd7UECTA>  .FO. I^IN^OStL «I-1) )  GC  TO  5 
GO  TO  A 

5  I G1=I 7 ( NE  CT A  421 
I  G2=I  7  ( NE  PT  3  4*) 

IG3=I7<NECTA  4t») 

I G4=I 7 ( SE  CT  A  4=5) 

IGPl=f 


IGP2=r- 
I GP3  =  t 
IGP4=L 

DO  6  J=1,L3«»LA 
K  =  NGPLA4J-1 


IFdGi.EO.I7  IK)  I 
IFIIG2.E0.I7  <K>) 
IFIIG3.E0.I7  IK)  ) 
IFd34.E0.I7  CK)) 
6  CONTINUE 

I  FdGPl  .20.  "  )  GO 
IF  (IGP2  .EO.  1 )  GO 
IF (IGP3.EO.  1 )  GO 


~GPi=NGPLA4tJ-l>4* 
■GP2=NGBIA4( J-l)44 
'GP3  =  NGPLA4(J-1)V‘- 
rG»'+=NG°IiA4<  J-l)4u 

TO  5 0 3 
TO  5*1 
TO  5  )  J 


IFdGPA.EO.O)  GO  TO  5 1 3 
7 (NV3LS45*I51-3) =1.3 
Z(NGVLS43*I32-7)  =1 . 0 
Z«NV3LS43*I33-3)  =1.0 
Z(NVGLS45*IG4-3)=-,0 

Z  (NVGMS45*  131-34  IT Gl-l)*N?OF)  =  1.0 
Z I NVGMS  45*132-74 (IG2-1) *N COF) =  1.C 


ZINVGMS 43  *133-74  CIG3-1 )  *N  fOF)  =1.0 
Z  CNVGMS 43 *134-34  f XG4-1) *N ?OF) =1.C 
X 1=Z  <  NBG? A4  r  G°l4l J 


Y  1=7 1 NGPA  *1 3P1*2  ) 


X2=7 1 NGPI  41  G°2*l ) 
Y  2=7 ( NGPA  41 5  P2*2> 
X  3=7  ( NGPA  4IGP341) 
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YSaZCNGP^ISP’*?) 

X4=7(KGfl*;3P4»l) 

X<*=7(NGPA+ISP4«-2> 

A  =  .?»  <<X3-X1)MY4-Y2>  +  (Y3-Y1)MX4-X2)  ) 
*ROW=KD 

oo  e  j=i,ro 

d  7  <NSftSEL*J-i>=  r 
7 (NSASEL+7-! >  =  7 

CALL  P*C<  ,  SA3"L  ,MC«SA> 

M  R0W=  L7  TT  TT 

no  n 

11  7 (NTTTTT  +  J-1)  =  1.  ' 

7  (NTTT7T  +  f  *tGl-7)=.  2^ 

7 <NTTTTT*f 1 rG2-3»=. 2* 

7 (NTTTTT  +  e  ’IG3-37-.  2^ 

Z(NTT77T*f 'IG4-3»=.2‘- 

CALL  PAC<«MMTTT*Tn  ,TTTTT  ,»C°B) 

7  CNNO  +2*1-1)  =  <Xi+X’*X3+XV  )/k, 

?.  (NNQ  +2*1  )=(  Yl  +  Y?+V3+Y'>)/'i. 

2 (NTEX-lf IC1)=X1 

7  (NTE  X-l»r,?)  =  X» 

7 <NTLX-1*IG?)=*Y 
7 (HT:y-i*iG4)=X4 
7  (N7£Y-l«-ISl>=*t 
y CNTEY-1*IG2>=Y2 
2  (N7£Y-1*I»?)='*,3 
7 (NT£Y-lfI34)=r4 
3  CONTINUE 

CALL  CLOyE(zCnfV 
CALL  CL  05  £  (SAP  £L  »1 ) 

CALL  CL  03  £  (TTTTT  *  t ) 

CALL  WFUFLCMCTA) 

CALL  WFTrPL(MCT71 
MC8A  ( 1 )  GL  5 
MCeA(2)«: 

MCBA(3)«LVGLS 
MCBA  ( 4 )  =2 
MCBA<5  >=t 
MCBA  (F  )=j 
NCBA(7)=; 

HCBB( 1) GMS 
MCBB(  2)  =* 

MC6B  ( 3 )  =NCOr 
MC9B  ( 4 )  =2 
MC8P  ( 5  )  =1 
MCBB(6)=J 
HCBB(7)=: 

NROWsLVGLS 

CALL  GCP£N(YGLS»7(BUP1),1  » 

CALL  PAOCCMNVSLSI  ,VSLS,MCPA) 

CALL  Cl  05  E(  if  6LS  »  1’ 

CALL  GOP£N(^G«S,  7C8«F2>  fl  j 
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o  o 


00  12  J=t,<3P 

CALL  PAC<C*eNVT •«<?«‘fw-l)*MtOFI  j VGKS»MC 39) 

12  CONTINUE 

CALL  CLOS  E  ( GHS ,  1' 

CALL  WPTr  FL ( MC3 A ) 

CALL  WPTTFU^rn) 

CALL  G0P-)’(»/TM,7fpuri)  ,i) 

CALL  GCPf  N<tfTR**<B'JFD  >1) 

NTE7"-  =  HSAFT. 

NT  6AR  =  f-  TEX ^♦N»° 

LTEXR  =  NG3 
LTEYR  =  NG3 
N  C=J 

C  LOOP  12  CALCULATE*  TH-  V T  *■  AND  VTR  MATRICIES 

00  13  J=1,)4GP 

I F  (7  (  NT  c<  ♦  J“l)  •  E°.  rtA  G )  GC  TO  12 
7<NTEXR+NC-t)s7(NTEX*J-l> 

7  (NTEYf-  +NC-n*7(MT;Y+J-l) 

NC=NC+1 

13  CONTINUE 
NROW=NC 
NCBA  (D^TI 
MC96(2)=. 

MC9M3)=NF04 
MCPA  (A  )=2 
MCnA(£)=l 
NCPAtfJs: 

MCpA  <  7  J  =„ 

MCBB  ( 1)  =/TP. 

NCB9  (  2  )  =.• 

MC9B<3)=NP1* 

MC9P(4)=2 
MC9A<E  )=l 

NCB8(7)=.: 

CALL  PAC<  (T(NT"XR1  ,VTN,MCBA) 

CALL  PACCITCNTEYR*  ,¥TR,  HC99) 

CALL  CL  05  F(  if  TH,  1 ) 

CALL  CLOSE (l/Tp»  1 ) 

CALL  HRTTPLCMC3A) 

CALL  WRIT FL (MCBBI 
MIL=1CICjC: 

1  =  1 

NCB  EQUALS  THE  N'WBER  OF  CHOROWISE  BOXES 
NS9  EQUALS  T«F  NMm°ER  3F  S°ANWISE  BOXES 
NCB=7  (NACPTfS) 

NSB=7  (NACPT*2> 

00  14  J=l,MSB 
DO  14  1=1, NCB 
IAl=MIL*l 
IA2=PU*L*1 
2A3=MILH+NCB 
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I  A4=ML+L  ♦NCB^l 
i  a  Di  = : 

I  A  P2  =  I 
I  A  P3  =  V 
I  A  PL  =  T 

DC  1>  KK=1,LG?LA 
<  =NG°L  A ♦<  K-l 

IFCIA  l.ED.I7  CO  >  7A=>i=N5FLAMKK-l)‘- 
I  F  (li  2  .  tD  »I7  CO  >  TA°2  =  N5*:LAM‘<K'-l)*-'( 

IF  (I  A  3.  ED.i7  (O  )  IA^7  =  'jr.F>LA  +  <KK-l)*iT 
IF(l'"*r1»T7  CO)  I  A=>'»  =  MGPLA*  CKK-1)*!» 

15  CONTINUE 

IF(IAPl.EO.l)  50  TO  F 1 
IF(lAP2.t0.3)  50  TO  5.M 
I F  (I  A  P3  «£(?•  ?  )  50  TO  5  10 
I F  (I  A  Ft  ,f.r.  0)  50  TO  5  1  0 
Xl=7(NGP„A«-rS°l*l> 

X2=7(NGFL  A«-IA»24-1> 

X  3=7  ( NGPL.  A  ^1 AP3 ♦  1’ 

X4=?(NGP_  A*I 
Y  1=7  <  NGPL  A  A°1  ♦  ?’ 

Y2=7(NGPLf  *1  S°?0' 

Y3=7(NGPlA*IAd3  +  ?' 

Y4=7(NGPlA*IAO'+*2') 

N  NI=NT  E  X 

THESE  ARE  THE  COORDINATES  OF  A  Ef  0  POINTS  AT  THE 
ONE  DUARTE7  CHJ?0  o0I.JT  0F  THE  POX 
7 (NN I ♦2*1-1)  =(Xl+X74(X2  +  X4-X3-Xl)/4. ) /2. 

Z  (NNI  +  2*L)  =  (  Y1*Y?*Y3>Y4)  /4. 

L  =L*1 

14  CONTINUE 

HC9A(i)=5TA< 

HCBA  (  2)  =1 
MCPA  (  3 )  =N  I 
HCBA  (4  )  =2 
HCBA  (5  )  *1 
MCBA  (6)  *i 
HCBA  <7 )  *3 

ngtak=nni*ni 

KT*t 

KX*« 

K  Y*0 

<0=1* 

D**l 

NSCR*  C  INI  (NTtT)  ♦  W*ND4M*  CNIt3)  ) 

THIS  CHECKS  OPEN  "OPE  TO  SEE  IF  ENOUGH  IS  AVAILABLE 
IFCNCORE.  LT.  INSDP*NGTA<))  GO  TO  Elf 

CALL  SSPL  IN  (  NI »  7  f  JNI)  ,N0,  T«NND)  , KX,  KY  ,  KO  ,  KT  ,  07,7  CNGTAO  ,NSCR» 
*  ISING) 

9UF1*LC0R  E“IBUe* 1 

BUFesBUFl-IBU17 

NCORE*eUF2-t 
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IP(N0Ci- E.  LT.  (NGTAr+NI'ND)  )  GO  TO  510 
CALL  GCPEMSTAK, '’"(9UF1)  ,  i  ) 

NROW=NO 
DO  Id  J=1,'0 

CALL  PAC<  <t(NC’tAKMJ-1)*NC),GTAK,MC9A) 

16  CONTINUc 

CALL  CLOSE(ST4<<r,  1) 

CALL  WATJTLfv^) 

CALL  GLPrt  (t/T*M,  ,l  1 

CALL  GCPEMVT31,*,CUF2»  ,1) 

NTFX  =  NGTAK 
LI  EX  =  2*  M 
NTEY  =  N7  EX  +LT  EX 

LOOP  17  Z ALCUL AT ' r  THE  VT NA  AND  VTPA  MATRICI ES 
DO  17  1=1,  SI 

7  (NTEX  +  I-1)  =  7(MTI*2*I-1> 

Z  CNTEY+I-l)s?(NMI*2*I) 

17  CONTINUE 
MC9A(1)=\/T  «A 
MC9A  (2)  =  ! 

HCnA  (  3 )  =vlFD H 
HCPA  (L )  =2 
NCSA  ( i  )  =1 
HC9AC6)=r 
rtC9A<7)=: 

MC09(l)=YTri 
HCBP(  2)  =. 

NC9B<3)=MF0* 

MC99(i»)=2 

HC99(5)=1 

MCB9<t)=S 

MC99(7)=: 

CALL  PACX  (T(NT“X)  t.VTHA,MCBA) 

CALL  PAC<(7(NT~Y),VT R1,MC«?9) 

CALL  CL03E(VTMt,i) 

CALL  CLOSE(VTRA,it 
CALL  WRTT RL  C  MCBA ) 

CALL  MRTTPL  CMC99) 

N  AIDMIT  =1 
DO  16  1=1, NI 
Z  INAIDKIT  ♦I”1)=1»0 

18  CONTINUE 

CALL  GO  PE  NCAIDHIT»Z(Rijfi)  ,d 

MC8A  <  1 )  =A  IOHIT 

HCBA  <  2  )  =3 

HCBA(3)=MI 

MCBA  (A )  =2 

MCBA  (5 )  =1 

HCBA  (6)  -? 

NCSA  (7  )  =i 
NROMsNI 

CALL  PAC<CMNAI0M1TJ,AIDMIT,MCBA) 


I 


CALL  CLOSE  (  AITHTT  ,  1) 

CALL  Wf..TrFL(Mrri ) 

CALL  GOPEMSA<5,-MpUPi) 

T  VPOUT=  3 
NSAS=i 
Nt'OW=Z*  HI 
NTL=HI  *  NRO* 

CO  1 9  1=1  ,  IT  L 
7<N  SAS  )=:■••© 

19  CONTINUE 
L  =  1 

00  2.'  1=1  ,NI 
Z<N$£S*L*I-?)=l,0 
7 <NSAS+L*I-l)=i. 0 
2f;  CONTINUE 
HC9A(1)=SAS 
MC9A  (2)=: 

NCPA<3)=MR0R 
MCP*  («< )  =2 
1C8A  ( 1  )  =3 
MC9A  ( 6 )  =  . 

MCeA(7)=. 

L  =1 

00  21  1=1, NT 

CALL  PA  C<  (T  (  NSA  9  ♦?  -1)  ,  SAS  ,MCr>A) 
L*L+NPCV» 

21  CONTINUE 

CALL  CLOS  E  { 5  AS  ,  1 ) 

CALL  WATT  PL  <  MC3  A  ) 

RETURN 
500  N=-fl 
FILE=C 

CALL  MES1GE(N,FILE,HA'1> 

END 
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